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Recently a nonuniversal character of the leading spatial behavior of the thermodynamic Casimir 
force has been reported [X. S. Chen and V. Dohm, Phys. Rev. E 66, 016102 (2002)]. We reconsider 
the arguments leading to this observation and show that there is no such leading nonuniversal term 
in systems with short-ranged interactions if one treats properly the effects generated by a sharp 
momentum cutoff in the Fourier transform of the interaction potential. We also conclude that lattice 
and continuum models then produce results in mutual agreement independent of the cutoff scheme, 
contrary to the aforementioned report. All results are consistent with the universal character of the 
Casimir force in systems with short-ranged interactions. The effects due to dispersion forces are 
discussed for systems with periodic or realistic boundary conditions. In contrast to systems with 
short-ranged interactions, for L/£ 2> 1 one observes leading finite-size contributions governed by 
power laws in L due to the subleading long-ranged character of the interaction, where L is the finite 
system size and £ is the correlation length. 
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I. INTRODUCTION 

According to our present understanding, the Casimir 
effect is a phenomenon common to all systems character- 
ized by fluctuating quantities on which external bound- 
ary conditions are imposed. 

The confinement of quantum mechanical vacuum fluc- 
tuations of the electromagnetic field causes long-ranged 
forces between two conducting uncharged plates which 
is known as the (quantum mechanical) Casimir effect 
0. & H The corresponding force between the plates 
is called the Casimir force. In this form the phenomenon 
was predicted in 1948 pj by the Dutch physicist Hendrick 
Casimir. 

The confinement of critical fluctuations of an order pa- 
rameter also induces long-ranged forces between the sys- 
tem boundaries 0,0,0- This is known as the statistical- 
mechanical (thermodynamic) Casimir effect. In this form 
the effect was discussed by Fisher and de Gennes al- 
ready in 1978. 

Casimir forces arise from the influence of one portion 
of a system, via fluctuations, on another portion some 
distance away. 

The best known example of a Casimir force is the van 
der Waals interaction between neutral molecules. In this 
case the correlations between fluctuations are mediated 
by photons, i.e., massless excitations of the electromag- 
netic field. When the system is a thermodynamic one, 
important examples of such massless excitations include 
Goldstone bosons and order parameter fluctuations at 
critical points. 

The quantum mechanical Casimir effect has been ex- 
perimentally verified with impressive experimental pre- 
cision (see also Refs.0 and 01 )■ One uses atomic 
force microscope techniques and measures the force be- 
tween a metallized sphere and a plate. Since it turns out 



to be very difficult to keep two plates parallel with the 
required accuracy, there is only one recent experiment 
[ill ] in which the original theoretical parallel plate set-up 
as studied by Casimir was used. In this experiment it has 
been found that the measured Casimir force agrees with 
the predicted one within 15% accuracy. It is interesting 
to note that at distances of the order of 10 nm between 
the plates the force produces a pressure of the order of 1 
atmosphere. Therefore, the Casimir effect is considered 
to be very important for the design of nanoscale devices 
(see, e.g., Ref.fT^] and references therein). 

In statistical mechanics the Casimir force is usually 
characterized by the excess free energy coming from the 
finite-size contributions to the free energy of the system. 
The parallel plate or film geometry turns out to be of 
great practical importance for experimental setups. 

A useful model for the investigation of generic finite- 
size effects is given by an 0(n)-symmetric spin system 
(n > 1) confined to a film geometry (L x oo 2 ) with pe- 
riodic boundary conditions r. Models of this sort serve 
as theoretical descriptions of magnets or fluids confined 
between two parallel plates of infinite area. The Casimir 
force per unit area in these systems is defined as 



^Casimir (-^5 



dL : 



where f£ x (T,L) is the excess free energy 

f^(T,L) = r(T,L)-Lf hulk (T). 



(1.1) 



(1.2) 



Here f T (T, L) is the full free energy per unit area (and 
per fcgT ) of such a system and /bulk is the bulk free 
energy density. 

According to the definition given by Eg. 1)1. If) the ther- 
modynamic Casimir force is a generalized force conjugate 
to the distance L between the boundaries of the system 
with the property F^ simiT (T, L) — > for L — > oo. We 
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are interested in the behavior of F^-- when L a, 
where a is a typical microscopic length scale. In this 
limit finite size scaling theory is applicable. The sign of 
the Casimir force is of particular interest. It is supposed 
that if the boundary conditions r are the same at both 
surfaces ^casimir wm be attractive 0, El E3 (strictly 
speaking, for an Ising-like system this should hold above 
the wetting transition temperature T w \13L Il4l Il7|). In 
the case of a fluid confined between identical walls this 
implies an attractive force between the walls for large sep- 
arations. When the boundary conditions differ between 
the confinin g su rfaces, the Casimir force is expected to 
be repulsive \13L VL3L [l6j . The current experimental situ- 
ation will be discussed later in the article. Here we only 
mention that these experiments are in qualitative and in 
some cases even in quantitative agreement. 

In this article we discuss the behavior of the thermo- 
dynamic Casimir force in systems with short-ranged and 
with subleading long-ranged (dispersion) forces, which 
are present, e.g., in real fluids. Both interactions lead 
to the same universality class, provided that the dimen- 
sionality d of the system and the symmetry of the or- 
dered state are the same. Despite this similarity we will 
see that, in comparison with systems with short-ranged 
forces, new important finite-size contributions exist in 
systems with dispersion forces. We shall also discuss 
proper boundary conditions for systems with subleading 
long-ranged interactions and ww shall reconsider several 
recent statements [l8| for the behavior of the finite-size 
free energy and the Casimir force in systems with short- 
ranged interactions and with dispersion forces. 

From the definition in Eas. l|l.l[) and (|1.2() it is clear 
that one needs to know the critical behavior of the free 
energy in a slab geometry in order to derive the behavior 
of the Casimir force. Based on numerous investigations 
it has turned out that the thermodynamic behavior of a 
system near a second order phase transition exhibits scale 
invariance and universality pH l2(ij | . In order to set the 
stage for our considerations we first recall certain bulk 
properties. 

A. Bulk systems 

Scale invariance and universality hold for the singular 
part of a thermodynamic function. For later reference 
we quote the decomposition into a regular and a singular 
part of the free energy / in units of ksT c and per unit 
volume of, e.g., an Ising ferromagnet: 

/(*. h) = f reg(t, h) + f sing (t, h) 

= freg(t, h) + \t\ 2 -'*A l F ± (A 2 h\t\~ A ) ,(1.3) 

where t = (T — T c )/T c ^ is the reduced temperature, h 
is the external magnetic field, A is the critical exponent 
associated with the magnetic field, a is the critical ex- 
ponent of the specific heat, Ai and A2 are nonuniversal 
(system dependent) metric factors, and F± are universal 
scaling functions. 



In the scaling limit the two-point correlation function 
in zero field, which is of particular interest in the present 
context, has the form [2l|, [23] 

G(r,t £ 0) = Dr-( d - 2 +^g ± (r/Z ± ) (1.4) 

with £± as the correlation length given by 

Z±(t,h = 0)=tf\t\-», t^O. (1.5) 

The nonuniversal constants D and £q can be related to 
5f, Ax, and A 2 via £+/£q = Qi, A Y = Q 2 ($)- d , and 

A 2 = Q3^D(^o) id+2 ^ )/2 with Qi> Q2, and Q 3 being 
universal, which leads to the hyperuniversality hypothe- 
sis in the form of two-scale factor universality [23| . In the 
form given above Eqs. (|1.4H and l|1.5fl are valid for Ising- 
like systems only. For 0(n) models, n > 2, one has in 
addition to take into account that £_(f < 0) = 00. The 
universal scaling function g + (x) decays exponentially for 
x > 1. The physical origin for the onset of scale invari- 
ance can be traced back to the divergence of the cor- 
relation length £±. Consensus has emerged that these 
statements hold in systems governed by short-range in- 
teraction potentials, i.e., decaying exponentially or being 
of a finite range. 

If the exchange interaction in an Ising model on a lat- 
tice in d dimensions decays algebraically, 

J (r)= 1 + (r / a)rf+CT > r^\r\>a>0, (1.6) 

where a is the lattice constant, the value of the decay ex- 
ponent a is crucial with respect to universality. For <r > 2 
the leading thermodynamic critical behavior is charac- 
terized by critical exponents and scaling functions for 
short-ranged interactions [24J. Mean-field theory holds 
for d > d c — 4 irrespective of the value of a. For a < 2 
the upp er critical dimension is reduced to d c (a) = 2a 
[24l |25| and the values of the critical exponents depend 
on a for a < d < d c (a) [2J, |2f| |23 . The crossover from 
short-ranged to long-ranged critical behavior occurs for 
(7 = 2 — r]sr(d), where rj sr {d) is the critical exponent for 
the short-ranged system (for a given fixed spatial dimen- 
sion d) [2^, I23, H3, El Ull • This crossover has recently 
been reexamined numerically in d — 2 in Ref. |33| . 

Fluids are governed by dispersion forces. In the sense 
of Ea. l|1.6[) dispersion (van der Waals) forces in d = 3 
dimensions are characterized by a = 3 in the nonre- 
tarded case and by a = 4 in the retarded case. Therefore, 
the leading thermodynamic critical behavior of a fluid is 
characterized by critical exponents and scaling functions 
for short-ranged interactions and the contributions due 
to the power law decay of the interaction potential lead 
to corrections to the asymptotic scaling behavior. Thus 
for a > 2 we refer to interaction potentials governed by 
Ea. (|l.tj|) as subleading long-ranged interactions. 

As an illustration of this case Fig. ^ displays schemati- 
cally the two-point correlation function G(r, t). The uni- 
versal decay of G(r, t) is governed by Eg . 1^1 .411 only within 
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FIG. 1: Schematic view of the density-density correlation 
function G(r, t) in a fluid governed by dispersion forces in 
d = 3. The behavior is shown on various length scales, where 
the tilted double slashes / / indicate breaks in scale. For mi- 
croscopic distances of the order of the particle diameter ao 
"packing" effects lead to oscillations decaying exponentially 
|34j . Beyond a crossover regime (not shown) the correlation 
function decays according to the power law given in Eq. ll.4ll 
as long as r <C £. Beyond another crossover (not shown) the 
further decay is exponential for r > 5 which finally crosses 
over to the interaction dominated regime for r > r* , where 
the ultimate decay of the correlation function follows the de- 
cay of the interaction potential ~ r -6 in the nonretarded 
regime and ~ r~ 7 in the retarded regime (not shown). The 
behavior r* ~ £ln£ of the crossover distance r* |35l l36l l37| 
illustrates that r* diverges more strongly than the correlation 
length f in the vicinity of T c , i.e., the critical regime, in which 
the universal properties hold, expands as T approaches T c . 



the critical regime. For distances r smaller than the lower 
limit of this regime generically nonuniversal microscopic 
effects govern the behavior of the correlation function. 
For distances larger than the upper limit of this regime 
the interaction potential itself governs the further de- 
cay of the correlation function. Note that r* /£ — > oo 
as T — > T c , i.e., the critical reg ime expands as the criti- 
cal point is approached j3f|, I36t l37j . The widening of the 
critical regime leads to the divergence of the compress- 
ibility /e(f) - r o °° drr 2 G(r : t) ~ |t|-f, for t 0. 

For short-ranged interactions one has to bear in mind 
that Eq. I|1.4|) is also only valid within a critical regime r < 
r*„. Exact results for the two-dimensional Ising model 
[33 and mean- field results (3^ suggest that r* r ~ £ 2 , 
whereas the spherical model yields the estimate r* r ~ £ 3 
[40|. For r 3> r* r the correlation function decays again 
exponentially, but it contains a nonuniversal prefactor 
|40j, i.e., the leading-order behavior becomes nonuniver- 
sal. This demonstrates that in the case of short-ranged 
interactions the width of the critical regime is much larger 
than for a corresponding system with subleading long- 
ranged interactions. 



B. Finite size scaling 

Finite-size scaling asserts that near the bulk critical 
temperature T c the influence of a finite sample size L 
on critical phenomena is governed by universal finite-size 
scaling functions which depend on the ratio L/£, so that 
the rounding of the thermodynamic singularities sets in 

for in ~ 0(1) HElElElElEi 

From the above discussion of the behavior of G(r, t) 
one expects that deviations from standard finite-size scal- 
ing behavior will be observed for L 3> r*, where r* is 
a crossover length with the property r* 3> £• In par- 
ticular one has r* ~ £ 2 for the Ising model or within 
mean-field theory, r* ~ £ 3 for the spherical model |4C| . 
and r* ~ £ log £ for subleading long-ranged interactions 

nam. 

This view has recently been challenged by Chen and 
Dohm [l8l | who purportedly report, for systems with 
short-ranged interactions, leading finite-size contribu- 
tions different from the ones expected from the above dis- 
cussion. This would invalidate the current understand- 
ing of finite-size scaling. In particular these results can 
lead to the expectation of a nonuniversal Casimir force 
at T c for a fluid between parallel plates at a distance L. 
If correct this would be of major theoretical 0, l48l |49| 
and experimental interest |5(J, |5ll |52] • Specifically, based 
on exact results for the 0(n) symmetric 4> 4 field theory 
in the large-n limit (mean spherical model), the authors 
of Ref. [I3 report the following result for the singular 
part f s (t, L, A) of the finite-size contribution f(t, L, A) to 
the free energy density of a system with periodic bound- 
ary conditions and purported short-ranged interactions 
in 2 < d < 4: 

/,(*, L, A) = L- 2 k d - 2 <S>{£- 1 h.- 1 )+L- d X ar (£/£). (1.7) 

The parameter A is an ultraviolet momentum cutoff and 
the function $ has the property $(0) > 0. For the 
Casimir force defined by 

F C a Si rmr(t, L, A) = - JL{i[/( t) L, A) - f(t, OO, A)]}, 

(1-8) 

Eq. (|1.7|) implies a leading nonuniversal (cutoff depen- 
dent) non-scaling term ~ L~ 2 in the behavior of the 
Casimir force, because the scaling function X sr (x) ~ 
exp(-x) when x > I 0, El IS El S3- There- 
fore Eas. l|1.7|l and (|1.8|) would also imply nonuniversal 
Casimir amplitudes 

A Casimir (d) = A d - 2 <£(0) (1.9) 

in d > 2 dimensions. 

In the following we shall show that the results reported 
in Ref. [18( can be traced back to using a peculiar model 
in which the interactions are neither short-ranged nor of 
the subleading long-ranged type so that the model does 
not relate to any physical realization. We find that if the 
periodicity and analyticity of the Fourier transform J(k) 
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of the interaction J(r) at the boundary of the Brillouin 
zone (in the case of a lattice model) and the analyticity 
of J(k) at the cutoff k = A (in the case of an off-lattice 
model) are preserved in the theoretical analysis, then the 
L~ 2 term in Eg. (|1.7l) vanishes identically. We also show 
that the presence or absence of this term does not depend 
on the range of the interactions. If the above require- 
ments for J(k) at the boundary of the Brillouin zone or 
at k = A are violated, a corresponding non-universal, 
non-scaling term of order L~ 2 will be observed in the 
finite-size behavior of any thermodynamic function. A 
discussion on the influence of the cutoff on the finite-size 
behavior of the susceptibility has already been presented 
in Ref. 0|; see also the "note added in proof" of Ref. 
|H^ . In Sections II and III we present a general and 
unified approach that is designed to avoid similar artifi- 
cial effects; this should be useful also in the context of 
quantum phase transitions and field theory. In Section 
III we also summarize the present state of knowledge on 
the finite-size behavior of systems with subleading long- 
range interactions, focusing on the expected behavior of 
the singular part of the free energy. One should distin- 
guish between the cases i) d + a < 6 and ii) d + a = 6, 
which contains the physically most important case of dis- 
persion forces in d = 3 with a = 3. In case ii) one expects 
additional logarithmic finite-size contributions. The ar- 
ticle closes in Sect. IV with a summary and concluding 
remarks, where we also discuss possible boundary con- 
ditions for systems with subleading long-ranged interac- 
tions. 



II. FINITE-SIZE BEHAVIOR OF THE FREE 
ENERGY DENSITY 

In this section we present our critique of the finite- 
size scaling analysis of the free energy and the Casimir 
force presented in Ref. . As pointed out already the 
statements of Ref. are based on exact results for the 
0(n) symmetric 4 field theory in the large-n limit (mean 
spherical model) with periodic boundary conditions. 



A. Analytical properties 

Before we turn to the finite-size analysis, we dis- 
cuss briefly the consequences of the assumptions for the 
Fourier transform of the interaction J(k) used in Ref. 
fl8| for the bulk properties of the model. For a system on 
a lattice J(k) = VJ r J(r) exp(ik.r), where the sum runs 
over the lattice sites. For an off-lattice system the sum 
has to be replaced by the corresponding integral. 

In Fig. |3 we compare J(k) for a short ranged (near- 
est neighbor) lattice model (see, c.f., Eq. l|2.I3|) ) with the 
standard k 2 spectrum in the infrared limit for short- 
ranged interactions in the spirit of Ref. 0] endowed 
with a cutoff A = n/a at the boundaries ±ir/a of the 
first Brillouin zone in one dimension. The restriction to 
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exact 




FIG. 2: Dispersion relation 1 — J(fc)/J(0) as function of k for 
nearest-neighbor interactions in one dimension (solid line) in 
comparison with a fc 2 -spectrum with a sharp cutoff at the Bril- 
louin zone boundary, for d = 1 dimension, as implemented in 
Ref. |lg(l (dashed line). The zone boundary is marked by the 
vertical dash-dotted line and a is the lattice constant. Note 
that the application of a sharp cutoff to a pure k 2 spectrum 
in the first Brillouin zone implies an artificial cusp-like non- 
analyticity of the dispersion relation at the zone boundaries. 
This is absent for genuine short-ranged interactions. 



one dimension only simplifies the notation and is not es- 
sential for the following arguments. ^From Fig. [2] it is 
obvious that the continuation of the fc 2 -spectrum, which 
is correct only in the infrared limit k — > 0, to the full 
Brillouin zone and its truncation at the zone boundaries 
introduces a cusp-like nonanalyticity into the spectrum. 
This nonanalyticity is artificial and is not a generic fea- 
ture of short-ranged interactions. As we shall demon- 
strate below, this artificial choice is the reason for the 
non-scaling and non-universal finite-size effects ~ L~ 2 
reported in Ref. 0|- I n fact, the importance of the 
properties of the dispersion relation at the Brillouin zone 
boundary for the critical finite-size behavior of the free 
energy in d = 2 dimensions has already been mentioned 
by Cardy (see Eq.(3.12) of Ref. [5f|) in the course of 
deriving the Casimir amplitude Ac a simir(2) — —irc/6 
in d — 2 for periodic boundary conditions, where c is 
the central charge of the model under consideration (e.g., 
c = 1/2 for the critical 2d Ising model with short ranged 
interacions). 

It is instructive to investigate the consequences of a 
nonanalytic spectrum of the kind displayed in Fig. [21 
in real space. According to Ref. [li| the corresponding 
correlation function reads for r 3> £ in d dimensions 



G(r,t) 



2A d - 2 (27rrA)-( d+1 )/ 
O(exp(-r/0) 



1)/2 sin[Ar-7r(d-l)/4] 



i + e- 2 A- 



(2.1) 



Thus the correlation function decays according to a 
power law with the decay exponent (d+l)/2 rather than 
exponentially. Furthermore, the correlation function os- 
cillates with a period set by the inverse of the cutoff A. 



5 



For separations r > r* (see Fig. Ea. (|2.1|) therefore 
implies that the interaction potential for a system with 
a truncated fc 2 -spectrum as shown in Fig. should be 
leading long-ranged rather than short-ranged or even sub- 
leading long-ranged in 2 < d < 4, but also containing 
both positive and negative contributions. Therefore, for 
all spatial dimensions of physical relevance, the model 
investigated in Ref. |l8j would appear to correspond to 
a model with competing leading long-ranged interactions 
in real space. For such a model the finite-size scaling as 
developed for systems with short-ranged or subleading 
long-ranged interactions is not expected to be applica- 
ble. 



B. Finite-size properties of the 0(n — > oo) model 

We substantiate our view by turning to a detailed anal- 
ysis of the finite-size behavior of the free energy and 
other thermodynamic functions, i.e., we provide an ac- 
count of the mathematical mechanism that produces the 
non-universal and non-scaling leading finite-size effects 
reported in Ref. |l8j . 

As a case study we quote the expression for the to- 
tal free energy density of a fully finite mean spherical 
model (the n — ► oo limit of an O(n) model) with nearest- 
neighbor interactions of strength J on a hypercubic lat- 
tice, which is given by |7j 



Pf(K, h\L, A) = ~ sup {--£- + — 



1 



t + ^2(1 — cos kj ) 

3=1 




K 

In 2dK 

2tt 



where K — (3J and ft is a properly normalized mag- 
netic field. The parameter r is related to the correlation 
length £ via r = £ -2 (4^ which relates r also to the re- 
duced temperature t. The relevant physical information 
of Ea. H2.2fl is contained in the sum over the wave vec- 
tors k = (fci, . . . , kd) in the first Brillouin zone of the 
simple cubic lattice. For general interaction potentials 
on general lattices this sum involves the dispersion re- 
lation w(k), where, e.g., w(k) = Ylj=i 2(1 — cosfej) for 
nearest-neighbor interactions on a simple cubic lattice 
(see Ea. H2.2fl and, c.f., Eg. H2.13f0 . In order to provide 
a general description of finite-size scaling of the free en- 
ergy, we will therefore consider the quantity |7J (see also 
Eq.(14) in Ref. 13) 



Ud,a{T, L, A) 



1 



2Li x 



x L d f-t 



ln[T + w(k)] (2.3) 



keB 



as function of the system size L, where L = (L 1; • • ■ , Ld) 
denotes the set of lengths that determine the geometry 
of the system, A = (Ai,--- , A^) is the set of cutoffs 



in reciprocal (i.e., k) space and r is a quantity propor- 
tional to the reduced temperature t, e.g., r = for 
the mean spherical model with short-range interactions. 
For lattice systems k belongs to the first Brillouin zone 
B and if the system is on a hypercubic lattice one has 
Aj = w/ai, where a, is the lattice spacing along the di- 
rection i, i — 1, • ■ • , d. For off-lattice systems the sum- 
mation is carried out over those values k€B that fulfill 
the requirements —A; < fc, < A,-, i = 1, ••• ,d. How- 
ever, for off-lattice systems, there is no obvious choice 
for the cutoff. One usually takes A ~ a -1 , where a is 
some fixed characteristic microscopic length of the sys- 
tem. For a fluid a can be taken to be the diameter do of a 
fluid particle (see also Fig. P). In Eq. I|2.3fl the subscript 
a characterizes the range of the interaction. The fluc- 
tuation spectrum u>(k) of the order parameter is given 
as a linear function of the Fourier transform J(k) of the 
interaction (see, c.f., Eas. 12.13fl and If2.14p l. 

The formal expression given by Eq. (f2.3f) . on which 
our specific considerations in this subsection are based, 
has widespread applications. The expression in Eq. I|2.3fi 
always appears as the one-loop contribution to the free 
energy in field-theoretic Ginzburg-Landau models |48| . 
The line of arguments presented here is therefore of gen- 
eral importance and is not limited to the specific model 
under consideration here. 

For periodic boundary conditions Liki/(2n) = Mi are 
integer numbers, i = 1, • • • , d, where Yli=i a Mi = N 
fixes the number of degrees of freedom in the system. The 
values of the components fcj of the vector k are given by 
ki = 2irmi/Li, with —Mi < < Mi — 1, i = 1, • • ■ , d. 

In order to analyze the sum in Ea. i|2.3fl we use the 
Poisson summation formula 

b oo „b -, 

E/H= E / dme i2 ^ n f(m)+-(f(a) + f(b)). 

m—a n— — oo a 

(2.4) 

After some algebra one obtains 

U dttT ( T , L, A) = U d , a (T, A) + AC/ d , CT (r, L, A), (2.5) 
where 

f-Ai /-Ad 



U dt<7 {r,A) = 



1 



dmi 



dnid ln[r + w(m)] 
(2.6) 

takes into account the contributions of the bulk system, 
while 



A C ^(r,L,A) = ^ r g£ i * Bl .../ 



drriri 



x exp I i^^njnijLj ln[r + o;(m)] (2-7) 

V j =1 / 

incorporates all contributions due to the finite size of the 
system. For further analysis of the Casimir effect the 
film geometry L x oo^ 1 is the most relevant one. It is 
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obtained from Ea. l|2.7|) in the limit L 2 — > 00, • • • , L d — > 
00, setting L\ = L. In order to simplify the notation we 
finally set Ai = A2 = ■ ■ • = = A so that 



m#0* 

x exp (inimiL) ln[r + w(m)], 



(2.8) 



which after two integrations by parts with respect to m\ 
can be rewritten as 



AU dA T 7 L,A) = AU d :i l(T,L,A) 



-L~ 



! ^FS4/-/ mi "7_ A A ' 



ni#0 

x exp (inim\L) d mi 



d mi uj(m) 



_t + aj(A,m 2 , - ■ ■ ,rn d ) 



(2-9) 



with 



Af7«(r,i,A) = L- 



dm 



6 (2^) Q j_ A 
d mi m(m)\ mi=A - d mi uj{m)\ mi= -A 
t + w(A,m 2 , • • • , m d ) 



e?m„ 



-A 



(2.10) 



The above expression is obtained by identity transforma- 
tions of the initial sum and is valid both for lattice and 
off-lattice systems. In the analysis below we show that it 
is the term AJjj^r, L, A) which produces the contribu- 
tions on which the statements of Ref. 01 are based. 

First, we evaluate this term for lattice systems. If the 
interactions J(r) > are such that they depend only 
on the distances between the particles, and these are 
the only interactions we are concerned with here, then 
J(k) = J(— k). Since there is no physical reason for sin- 
gularities anywhere except at k = the derivatives of 
J(k) with respect to k should exist at least for all k 7^ 
and therefore <9k</(k) = — c\J(— k). This holds for lat- 
tice and off-lattice systems. For lattice systems J(k) is a 
periodic function with the property 



J(k + 2A iei ) = J(k), 



(2.11) 



where is a unit vector in reciprocal space. This implies 
that <9k>/(k) = at the borders of the Brillouin zone and 
therefore 



AE#>(t,L,A) = 0. 



(2.12) 



Note that the above result does not depend on the range 
of the interaction - it is true for short-ranged, subleading 
long-ranged, as well as for leading long-ranged interac- 
tions. As an illustration of the above general arguments 
we recall that the exact Fourier transform of the near- 
est neighbor interaction on a <i-dimcnsional hypercubic 
lattice reads 



J(k) = cos kj = J (2d- w(k)) 



(2.13) 



where 



w(k) 



d 

E 



2(1 — cosfcj). 



(2.14) 



The general properties of J(k) that we have discussed 
above can be easily verified from Eqs. (|2.13l) and l|2.14|l . 

We now reconsider the quantity A[/^(t, L,A) if the 
exact spectrum is replaced by its asymptotic form, valid 
in the infrared limit k — > 0, for all k e B. This is a 
very common procedure in the theory of critical phe- 
nomena based on the general idea that only long wave- 
length (small k) contributions are important for the crit- 
ical properties of the system. This leads to w(m) = m 2 , 
which is the spectrum used in Ref. 0] . One immediately 
obtains 



A[/«(t,L,A) 



A 



A 2tt 2 1 

1 



dm? 



drrir, 



A" 



- A 2 + n 
1 1 



\d-2 



L 2 6(2tt, .; i 
1 



m 2 

1 ,1 
dmi ■ ■ ■ I drrid 



1 + t/A 2 + m\ H \~rnl' 



(2.15) 



We recall that in the spherical limit of the 0(n) model 
one has t = £~ 2 . Ea. 12.15fl exactly reproduces the 
nonuniversal leading but non-scaling finite-size contribu- 
tion to the free energy as reported in Eq.(16) of Ref.[l8| 
for the corresponding field-theoretical model. Similar 
contributions exist also for subleading long-ranged in- 
teractions to which we turn in Sect. III. According to 
the above considerations such nonuniversal (cutoff de- 
pendent) contributions of the order of L~ 2 will always 
appear if 



9J(k) 



dki 



dJ(k) 



k 1= A 



dkx 



(2.16) 



fci = -A 



Note that only the properties of the Fourier transform 
J(k) of the interaction at the boundary of the set B of 
allowed k values is important here. For a field-theoretic 
model these are definitely a matter of definition. For 
lattice models these properties follow automatically. Ap- 
proximating the spectrum w(k) by its infrared asymp- 
totic behavior leads to an artificial cusp-like singularity 
at the border of the Brillouin zone as is illustrated in 
Fig. [3 A corresponding approximation for Eq.(3.11) in 
Ref. 55] would lead to an incorrect prediction of the 
critical finite-size contribution to the free energy (see 
Eq.(3.12) in Ref. |55| ) leading to a vanishing Casimir am- 
plitude. 

Before we consider how to modify the definition of the 
continuous field-theoretic model as to avoid a nonzero 
AJ7^^(r, L, A), we make some general remarks. First, the 
considerations presented above can easily be extended to 
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any geometry of the type L 



d-d' 



X oo 



<d' <d. Sec- 



ond, in the above discussion we did not specify the type 
of the interactions - short-ranged, leading long-ranged, 
or subleading long-ranged. This implies that the L~ 2 
corrections in question exist for any type of interaction 
for periodic boundary conditions, provided Ea. <2.16[) is 
valid. Furthermore, further integrations by parts yield 
additional contributions of the order L~ 4 , L~ 6 , etc. In 
2 < d < 4 dimensions only the term L~ 2 is important, 
but in d > 4 spurious finite-size terms of the orders L~ 2 
and L~ 4 will be generated. Finally, we note that the 
precise form of the integrand in Eq. I|2.8[) was not used 
in the above analysis. This implies that spurious L~ 2 
corrections also occur in other quantities such as the sus- 
ceptibility _JE3, 13] (see also the discussion in Appendix 
G of Ref.|46|]), the specific heat |56|, etc. The influence 
of different cutoff types (sharp or smooth) and of the 
truncation of the expansion of the Fourier transform of 
the interaction on the finite-size behavior of the suscep- 
tibility has been considered in detail in Ref.[46j. The 
only difference with respect to the free energy is that the 
susceptibility diverges as L 1 ^ at T c leaving any L~ 2 con- 
tribution as a small correction, whereas the singular part 
of the free energy behaves as L~ d and therefore the cutoff 
dependent term of the order L~ 2 becomes dominant in 
the critical region for d > 2. 

Since there is no physical reason for the introduction 
of a sharp cutoff in the /c-space representation of a field- 
theoretic model, one option to avoid artificial L~ 2 con- 
tributions is to implement of a smooth cutoff. Various 
forms of smooth cutoffs are possible 0, 0] . For exam- 



ple in Ref . |54| a modified continuum Ginzburg-Landau 
Hamiltonian has been considered (see also Ref. 



:g-Lan> 
,57]): 



H = 



1 2 

2 r0tp 



2^ 



u (tp 2 



(2.17) 

The last term in Ea. l|2.17fl introduces the smooth cutoff 
which is parameterized by a wave number A. The finite- 
size effects of thermodynamic quantities differ substan- 
tially for the above Hamiltonian and the standard one 
with a sharp cutoff: In the framework corresponding to 
Eq. (|2.17|l thermodynamic quantities approach their bulk 
value exponentially as function of L for T T c fixed, 
whereas for the standard Ginzburg - Landau Hamiltonian 
with a sharp cutoff the bulk limit is reached according to 
the power law ~ L~ 2 for any temperature. In particu- 
lar, this has been observed for the susceptibility [53, the 
specific heat [5^, and the free energy 01 • The effect of 
a smooth cutoff prodecure on the finite-size behavior of 
the susceptibility was also discussed in detail in Ref. ^(| , 
where the aforementioned exponential finite-size behav- 
ior was recovered. As already noticed in Refs. |4jj, 
and [H^ the presence of a sharp cutoff is mandatory for 
the occurrence of the aforementioned nonuniversal L~ 2 
contributions to finite-size scaling. For general d > 2 it 
was also realized that a close relationship exists between 
a non-exponential large-distance behavior of the bulk cor- 



relation function (generated by the sharp cutoff) and the 
power-law finite-size behavior of both the susceptibility 
above T c |4jj and the singular part of the free energy 
|l8j. Here we have demonstrated that all such finite-size 
effects arise from a single origin and are unphysical math- 
ematical artifacts due to the imposed singularity of J(k) 
at the boundary of the allowed fc-values. This, in turn, 
generates long-ranged correlations in real space. In order 
to eliminate spurious finite-size contributions we propose 
the replacement (see Eq. 112.81) ) 

Af/ rf)CT (r, L, A) ACV(r, L, A) - AU^l (r, L, A) 



wk 3 ^ L dmi ---C dmd 

x exp (in\m\L) ln[r + w(m)] 

7r 2 1 f A f A 
-L~ 2 — — ; j I dm,2 ■ ■ ■ I dm 



(2.18) 



6 {2nYJ_ A j_ A 
^ d mi u>(m)\ mi=A - d mi oj(m)\ mi= - A 
t + w(A,m 2 , ■ ■ ■ ,m d ) 

and corresponding replacements generated by derivatives 
of Eq. (|2.18l) with respect to the parameter r in the def- 
inition of each model system regardless of the imple- 
mentation of a sharp cutoff. Within such a scheme the 
well established methods for field-theoretic calculations 
in the presence of sharp cutoff are preserved. For lat- 
tice system this is an identity transformation, because 
A{7j2(r,i,A) = as expounded above. Note that the 
replacements given by Eci. (|2.18J) and its derivatives with 
respect to r do not interfere with the treatment of bulk 
systems, since A[7^(r, L = 00, A) = 0. They only be- 
come important for studies of the finite-size scaling be- 
havior of systems endowed with a sharp cutoff. These 
replacements remove the artificial cutoff dependent finite- 
size contributions to thermodynamic quantities in the 
spherical limit n — > 00 of 0(n) models and to one-loop 
order for 0{n) models with finite n. 



III. SYSTEMS WITH SUBLEADING 
LONG-RANGED INTERACTIONS 

First, we briefly recall the finite-size behavior of sys- 
tems with subleading long-ranged interactions. 

In Ref. [43 and in Refs. |S7j and |5J| it was shown 
that the susceptibility of a finite system with dispersion 
interactions, which decay as r~ d ~ a for large distances, 
can be written for 2 < d < 4, 2 < a < 4, and d + a < 6 
in the form 

X (t,L) - L'^X^L/^bL 2 -^) (3.1) 
~ [X£ (L/0 + btf-o-iX 1 ; (£/£)] , 

where 

X s x r (x -» +00) ~ X^'+x-il" + 0(exp(-const. x)). 

(3.2) 
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For the long-ranged part one has 

X lr (x -> +00) ~ XL r 1 x-^/ u+<T + X lr 



-2-i/u-d 



(3.3) 



The amplitude b is a nonuniversal parameter that can 
be determined from the Fourier transform of the in- 
teraction. The first term of the asymptotic behavior 
of X£(x) yields the bulk corrections to scaling as pre- 
dicted by Kayser and Raveche j3f|, while the second 
term yields the leading finite size correction to the sus- 
ceptibility for L/£ 3> 1. This second term leads to 
X (t,L) - x(t,oo) ~ t~ dv -^ L-^ d+a \ LIZ > 1, i.e., the 
finite-size corrections to the bulk behavior are governed 
by a power law rather than an exponential function of 
the system size L. Finally, we note that for the physi- 
cally most important case d + a = 6 (e.g., d = a = 3 for 
non-retarded van der Waals forces in d = 3) one finds 
additional logarithmic corrections in Eq.fl^HJ [I3- EH- 
which can be incorporated by the replacement X!£(x) 

-> X%> 1 (x)hiL + X%> 2 (x) M- 

The behavior of the susceptibility outlined above is 
consistent with the behavior of the bulk pair correlation 
function in systems with dispersion forces |37| 



G(r,t) 



[9±(r/0+r 



— (cr— 2+77) Ir 



g%{r/Z)\. (3.4) 



The modified Fisher-Privman |60| finite-size scaling hy- 
pothesis for the free energy density in such systems can 
be cast into the form 

f s (t,L) = L- d X(L/tbL 2 -°-<>) (3.5) 
~ L- d [X sr (L/Z) + bL 2 - a - r 'X lr (L/Z)], 

where X sr (x — > +00) ~ X sr ' + x d + 0(exp(— const, x)) is 
the short-ranged contribution. For the long-ranged con- 
tribution one expects X lr (x) ~ X[ r x d+a+, i- 2 + X l 2 r x^~ 2 . 
The first term in the asymptotic behavior of X lr yields 
a new bulk correction to scaling that is due to the sub- 
leading part of the interaction (analogous to the corre- 
sponding terms predicted by Kayser and Raveche for the 
susceptibility |35j and observed in spherical model cal- 
culations |4q[). Its temperature dependence for T > T c 
is given by i*H-(<r+r)-2)i/_ For > 1 the second term 
leads to a finite-size contribution of the form L~( d+a \ 
As for the finite-size scaling behavior of the susceptibil- 
ity, additional logarithmic corrections have to be added 
to X lr for d + a = 6. Finally, we note that there may 
also be a third constant contribution to the asymptotic 
behavior of X lr (x) for x ^> 1, which does not appear in 
the susceptibility, and which leads to a L~ A In L finite-size 
contribution to the free energy for d = a = 3. 

Eq. 1)3.5(1 has been derived in Ref. 18j for the case r\ = 
where the main focus is set on the discussion of finite- 
size contributions. Furthermore, the results reported in 
Ref. |18| apparently apply only for d + a < 6, because 
no logarithmic corrections were found. Thus a complete 
verification of Eq. (|3.5|) is still missing. We now turn 
to the investigation of some of the consequences of the 
assumptions used in the model Ref. |Ts[ . 



We suppose that the interaction potential J(r) is of 
the dispersion type as defined by Eq. 1(1.61) . The Fourier 
transform of such an interaction is 



J(k) 



J(0) 
J(0) 



{l~v 2 k 2 + ; 



v 4 k 4 + 0(k 6 )) 



(3.6) 



where k = |k|, 4 > a > 2 and J(0), v 2 , v ai and V4 
are nonuniversal positive constants. The constant J(0) 
is the ground state energy of the system and w(k) ~ 
k 2 -bk a +ck 4 +0(k 6 ), where K = [3v 2 J{0), b = v tr /v 2 > 
and c = V4/V2 > are nonuniversal constants. Since 
J(r) > and thus J(0) > J(k) for k ^ 0, the values 
of b and c are such, that there are no real roots of the 
equation 1 — bk a ~ 2 + ck 2 = with respect to k. 

The free energy of an 0(n) model with an interaction 
described by Ea. (|1.6|) . in the limit n — > 00, is given by 
the expression 



Pf(K, h\L, A) = I sup I 

2 r>0 I i^T 

x H T + <*>(k)] - Kt 



1 



L\ x • • • x Ld 



1 



, K K 
In 

2tt v 2 



.(3.7) 



In the presence a sharp cutoff in fc-space this leads to 



AC#>(r,L > A) 



A d 1 



L 2 6 (2?r) 



d-2 



dm? ■ ■ 



drrid (3.8) 



1 - i6crA CT - 2 (l + e 2 )^ 2 - 1 + 2c 2 A 2 (l + e 2 ) 
t + A 2 (l + e 2 ) - 6A CT (1 + 6» 2 )°-/ 2 + 2c 2 A 4 (l + 6 2 ) 2 ' 



where 2 = + • 



■ + m 2 . 



This term is missing in Eq.(9) 



in Ref. |18j but it is manifestly present in the case of a 
sharp cutoff. We therefore conclude that Eq.(9) of Ref. 
|18| is only correct within the smooth cutoff procedure, 
but not within the sharp cutoff one. As already explained 
above, Eq. (9) of Ref. [3 coincides with our Eq. 1)3.5(1 
for systems with r\ = and d + a < 6 for periodic bound- 
ary conditions. In Ref.^| it is supposed to be valid also 
for systems with Dirichlet boundary conditions. How- 
ever, Dirichlct boundary conditions are inconsistent with 
the long-ranged nature of the dispersion forces (sublead- 
ing long-ranged interaction), because the "missing neigh- 
bors" of the ordering degrees of freedom at a surface of 
such a system by the nature of long-ranged interactions 
generate a long-ranged surface field. We therefore con- 
clude that the consideration of systems with subleading 
long-range interactions combined with Dirichlet bound- 
ary conditions as proposed in Ref.|l8j is of no physical 
relevance. In the next section we summarize our findings 
and also comment on the proper boundary conditions 
and the expected finite-size behavior of systems with dis- 
persion forces and real boundaries. We will present ar- 
guments as to why we expect this to differ significantly 
from Eq. ((3~5|) . 
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IV. SUMMARY AND CONCLUDING 
REMARKS 

It has been shown in Sects. II and III that nonanalyt- 
icities in the dispersion relation w(k) at the momentum 
cutoff lead to a bulk model which has leading and com- 
peting long-ranged interactions in real space. Thus even 
in the bulk it has peculiar properties such as the two- 
point correlation function given by Ea. (|2.1|l . For such 
a model finite-size scaling developed for systems with 
short-ranged or subleading long-ranged interactions does 
not apply from the outset. We are not aware of any 
physical system governed by such a type of interaction. 

In order to investigate the critical behavior of a model 
system by means of an effective Ginzburg-Landau Hamil- 
tonian in k-space, one usually expands w(k) in the (in- 
frared) limit k — > 0, keeping only the leading term(s). 
Divergent momentum integrals can be regularized by the 
introduction of a cutoff which is motivated by the pres- 
ence of a Brillouin zone (lattice models) or finite particle 
sizes (continuum models). If a sharp cutoff in momen- 
tum space is applied as a strict feature of an otherwise 
approximate (expanded) dispersion relation, the ensuing 
nonanalyticity of w(k) at the cutoff generates long-ranged 
correlations and finite-size effects of the order L~ 2 . These 
effects do not occur in actual systems with short-ranged 
or subleading long-ranged interactions and therefore they 
are artificial. In particular, the nonuniversal long-ranged 
Casimir forces reported in Ref. 01 have no physical rel- 
evance for systems with short - ranged forces. 

It is instructive to consider the finite-size behavior of 
the free energy in systems with subleading long-ranged 
interactions. As is well known, the free energy decom- 
poses into a sum of a regular and a singular part. In 
Section III we have discussed the finite-size behavior of 
the singular part for the case of periodic boundary con- 
ditions in a film geometry. However, the regular part 
is also important and it has experimental consequences 
for the Casimir force. For periodic boundary conditions 
one expects the regular part to be of the order L~ 4 in 
the vicinity of T c . Much more interesting is the case of 
a system with real boundaries. This raises the question 
of proper boundary conditions for such systems. They 
cannot be boundary conditions of the Dirichlet or Neu- 
mann type, because the latter are incompatible with the 
long-ranged nature of the interactions. Instead, one finds 
that long-ranged interactions generate long-ranged sur- 
face fields which decay according to a power law away 
from the surface into the bulk of the system. The direct 
(Hamaker) interaction of the surfaces then generates a 
L~ a contribution to the regular part of the free energy if 
the free energy is measured per unit volume (or a L~ a+1 
contribution if the free energy is measured per unit area). 
This is well known from studies of, e.g., wetting phenom- 
ena ^3 • The contribution to the regular part of the free 
energy due to the action of the surface fields on the or- 
dering degrees of freedom is also of the order of L~ a . The 
available renormalization group arguments suggest that 



the contribution to the sing ular part should be of the or- 
der £-(d-2-H)/2- ff |6ll|6||. This leads to the following 
hypothesis for the singular part of the free energy: 

f s (t,L) = L- d X(L/^h 1 L Al / u ,bL 2 -' T -'\h s L ( - d+2 - r i^ 2 - a ), 

(4.1) 

where h s is a nonuniversal metric factor characterizing 
the long-ranged behavior of the surface fields, while h\ is 
the corresponding factor at the surface boundaries. Here 
Ai is the critical surface gap exponent of the correspond- 
ing surface universality class. Note that, since d > 2, the 
term proportional to h s will give contributions that are 
larger than that proportional to 6, i.e., the contributions 
of the surface fields are important and cannot be ignored. 
In particular, this is important for the quantitave anal- 
ysis of wetting experiments with critical binary liquid 
mixtures |50| . The interpretation of these experimental 
data is still unresolved. We also point out, that nei- 
ther Ea. l|3.5fl nor Eq. (|4.1|l applies to confined superfluid 
4 He or 3 He- 4 He mixtures |48j which have been under in- 
vestigation in recent experiments 5l||52j. For superfluid 
4 He the long-ranged finite-size contributions to the free 
energy originate from two distinct sources: (i) a regular 
contribution from the dispersion forces which couple to 
the fluid density and which are unrelated to the super- 
fluid order parameter, (ii) a singular contribution from 
the superfluid order parameter which is generically short- 
ranged in nature and does not relate to the presence of 
dispersion forces ^||- The situation is more complicated 
in the case of 3 He- 4 He mixtures, where the superfluid 
transition temperature depends on the 3 He concentra- 
tion, which eventually causes the superfluid transition to 
become first order beyond the tricritical 3 He concentra- 
tion. The 3 He concentration does respond to dispersion 
forces and, by virtue of long-ranged surface fields, long- 
ranged 3 He concentration perturbations may emerge in 
a confined 3 He- He mixture. Through the dependence 
of the superfluid transition temperature T\ on the He 
concentration a long-ranged variation of the local value 
of T\ will ensue which in turn imposes a correspond- 
ing variation of the superfluid density in thermal equi- 
librium. From the experiment in Ref.[52j| there is robust 
evidence that Dirichlet boundary conditions do not apply 
for the superfluid order parameter of 3 He- 4 He in the tri- 
critical regime, because contrary to the case of pure 4 He 
pHj a thickening of the wetting layer has been observed 
caused by a repulsive tricritical Casimir force. This ob- 
servation rules out pure Dirichlet boundary conditions in 
this case, because these can only account for attractive 

Casimir forces as observed in pure 4 He mmi 

Further- 
more, subdominant long-ranged interactions may play a 
significant role in a finite-size scaling analysis of 3 He- 4 He 
mixtures in the vicinity of the bulk superfluid transition. 
However, the model Hamiltonian of the system then has 
to accomodate a second 'noncritical' field (the 3 He con- 
centration) apart from the superfluid order parameter in 
order to include dispersion forces and long-ranged surface 
fields in the physically correct way. The construction of 
a proper model Hamiltonian, which is a generalization of 
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the standard Ginzburg - Landau Hamiltonian considered 
here, is beyond the scope of this work. 

In conclusion, we remark that we still lack a complete 
theoretical description of the Casimir effect in systems 
with sublcading long-ranged interactions. Such a descrip- 
tion must contain both the influence of surface fields and 
the long-ranged nature of the interaction potential. Both 
are expected to generate important contributions to the 
critical behavior pertaining to the universality class of 
systems with short-ranged interaction potentials. Very 
close to T c the finite-size behavior should turn out to be 
that of short-ranged models, but additional finite-size ef- 
fects are expected to become dominant for \t\L}l v ^> 1. 
For systems with scalar order parameter these expecta- 
tions hold both above and below the bulk critical tem- 
perature. For 0(n) models one expects that Goldstone 



modes will dominate the finite-size behavior of the sys- 
tems below T c . 
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